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This papem considers the problem of estimating the time auto-correlation function for a quantity 
that is defined in configuration space, given a knowledge of the mean-square displacement as function 
of time in configuration space. The problem is particularly relevant for viscous liquids, where the 
interesting time-scales are often beyond those reachable by computer simulation. An approximate 
formula is derived which reduces the calculation of the time auto-correlation function to a "double 
canonical" average. In this approximation, in the case of Langevin dynamics the mean-square 
displacement itself may be evaluated from the "double partition function" . The scheme developed 
^v ■ is illustrated by computer simulations of a simple one-dimensional system at different temperatures, 

showing good agreement between the exact time auto-correlation functions and those found by the 
. ^ , approximation. 

(D 

QQ The calculation of a time autja-cprrelation functiorJj □ is a straightforward matter in any computer simulation tracing 

T-H the time evolution of a systenJj 13. However, computer simulations are not feasible today on time scales longer than 
microseconds. These time scales are relevant, e.g., for viscous liquids approaching the glass transition. Therefore, one 

^"^ ' cannot by simulation calculate_a number of experimentally easily accessible quantitiesinJijghly viscous liquids like the 

Jft frequency-dependent viscosityl3|-|bulk modulusQ, dielectric constanlll3, or specific healliJO (quantities that all via the 

^»^ fluctuation-dissipation theoremoQ are given as Laplace transforms of an equilibrium time auto-correlation function. 

*vj In this situation one would like to have an approximate theory at hand. Focussing only on time auto-correlation 

fvq . functions of quantities A{X) that are defined in configuration space, X = (Xi, ..., X„), an approximation is proposed 

T— I ' below, based on an ansatz for the joint probability of initial point X sXt — Q and final point X' after time i, denoted 

C ■ , by P{X,X';t). This joint probability directly determines the time auto-correlation function via 
0^ . 

^! {A{fi)A{t)) = f dXdX'A{X)A{X')P{X,X':t). (1) 

B- . ...... 

I ' If Z is the configurational partition function, given in terms of /3 = \/{kBT) and the potential energy U{X) as 

§: ZiP) - jdXe-^^i^\ (2) 

o ■ -^ 

^ , and G{X -^ X';t) is the Green's function, the exact expression for the joint probability is 

•I— I , 

H' P{X,X'-t) = ^-^—G{X^X';t). (3) 

However, Eq. (^ is not very useful unless the Green's function is known. Note that the principle of detailed balance 
- expressing the time-reversal invariance of the equations of motion - implies that P{X,X';t) — P{X',X;t), a 
requirement that any approximation should also satisfy. 

The exact method for calculating the time auto-correlation function {A{0)A{t)) is illustrated in Fig. 1 showing 
the path in configurational space. At a number of times ii, ...,tp one computes the product A{X{tj))A{X{tj + i)), 
and {A{0)A{t)) is the average of this quantity as p — > oo. Assuming here and henceforth that the X^'s are simple 
rectangular coordinates and that {A) = 0, one always expects {A{0)A{t)) ^ as i ^ oo. This loss of correlation after 
long time comes about because the point in configuration space at time t, X' , for large t is far away from the initial 
point X. A measure of the distance travelled in time t is provided by the mean-square displacement, (AA'^(i)). This 
quantity, of course, is defined by 

(AX^t)) = {(X X'f) ^ J2 ((^^ - X^)') ■ (4) 

4=1 

Assuming that the mean-square displacement itself is a known function of time, the idea is now to estimate 
{A{0)A{t)) via the "spatial" auto-correlation of A in configuration space evaluated at a distance equal to y^{AX'^{t)). 



Before proceeding, we briefly discuss the physics of this way of thinking about the time auto-correlation function. A 
simple case is when the mean-square displacement is proportional to time (for large times this is, of course, always 
the case). In this case, if the "spatial" correlation of A has a Gaussian distance decay, the time auto-correlation 
function is a simple exponential, corresponding to Debye relaxation. If, however, the spatial correlation of A has 
an exponential decay, the time auto-correlation function is a stretched exponential with exponent 1/2. The latter 
case gives a reasonable fit to many experiments on viscous liquidsE3. The above picture of decomposing the time- 
autocorrelation function into A) a "geometric" correlation and B) the distance travelled in a given time, is in harmony 
with another well-known property of viscous liquids. In these systems all linear relaxation functions determine roughly 
the same average relaxation time, a time which increases dramatically upon cooling. In the above "geometric" picture 
this is simply a consequence of the motion slowing down in configuration space, whereas "spatial" correlations in 
many cases probably only change little upon cooling. 

We now turn to the problem of estimating the joint probability P{X, X'; t). In the thermodynamic limit (?i — > oo) 
the relative fluctuations in the mean-square displacement go to zero, and therefore the distance between X — X{tj) 
and X' = X{tj + 1) is "precisely" equal to •\/(AAr2(t)). Similarly, the relative fluctuations in potential energy go to 
zero, so the potential energy of both points X and X' is "equal to" ([/} = — ^ L"j^ ■ Now, the basic idea is an ansatz for 
P{X,X']t) that assumes equal probability for all pairs of points with correct distance and correct potential energy. 
Thus, 

P(X, X'- 1) (X 5[{X - X'f - (AX2(t))] 6[U{X) - {U)] 5[U{X') - ([/)] . (5) 

In the thermodynamic limit there is " equipartition" between U{X) and U{X'), and the last two delta functions may 
be replaced by a single delta function, leading to 

P{X, X'; t) ex 5[{X - X'f - (AX2(t))] 5[U{X) + U{X') - 2 (C/)] . (6) 

The next step is to convert Eq. (|6|) to a "canonical" form, which is computationally much more convenient than the 
present " microcanonical" form. This is done by replacing the first delta function by cxp[— a(X — X'Y] where a is a 
Lagrangian multiplier adjusted to give the correct mean-square displacement. Similarly, the second delta function is 
replaced by exp[— &[[/(X) -I- [/(X')]] where h is adjusted to ensure that the average of U{X) + U{X') is 2 ([/). If the 
following "double partition function" 

D{a,b) = /dXdX'e-'^(^-^')'-''[^(^)+^(^')l (7) 



is introduced, the final expression for approximately calculating the time auto-correlation function is 

{A{0)A{t)) = /^54^A(X)A(X')e-"(^-^')'-''[^(^)+^(^')]. (8) 

J D{a,b) 

In the thermodynamic limit Eq. (pi) is equivalent to calculating the average of Eq. ([l|) by using Eq. (0) . 

The two parameters a and b are determined in the following way. First, b = 6(a) is found from the condition that 
the average joint potential energy of initial and final point is 2 ([/). Thus, b = b{a) is determined from the condition 
that this average is independent of a: 

^^^'^^ :0. (9) 



da db 



Since j^ ^ da + ^db [with standard abbreviation for partial derivatives], the expansion of Eq. (O) leads to the 
following first order differential equation for b = b(a): 

db_ ^ dgDdbD - DdabD 

da DdlD~{dbDY ■ ^ ' 

Once the function 6(a) has been determined, a = a{t) is found from requiring that the mean-square displacement 
calculated from D{a,b{a)) is correct: 

-^.(AA-(,)>^ (11) 

The short and long time limits are determined as follows. For a{t) one clearly has 



ait = 0) = oo and a(oo) = 0. (12) 

In the limit of large times X and X' are far apart and U{X) is uncorrelated with U{X'). In this limit h = (3: 

6(a = 0) = fi. (13) 

In the short time limit the points X and X' are close. Thus, P{X,X']t) oc 6[X — X'] exp[— 2feL/(X)] for t — > and 
Eq. (1) yields 



lim(A(0)A(t)) = 



For this to give the correct canonical average one must have b — /3/2, i.e., 

5(a = oo) == - . 



(14) 



(15) 



As we will show now, the short time behavior of a{t) may be derived directly from the equations of motion. In 
the case of Newtonian dynamics, the Green's function at short times is easily found from the integrated equations of 
motion, where the momentum is distributed according to a Gaussian (for simplicity all particles are assumed to have 
equal mass, m) 



G{X -^ X';t) ex exp 
To first order in t this becomes 



.JL 

2m 



X:(y(x;-xo + ia,t/t 



where 



G{X -^ X';t) ex exp 



j3m 



w^2 /3r 



-a{t){X-X'f-'-;^[U{X')-U{X)] 



a{t) = — 2" [Newtonian dynamics, i — > 0] 



(16) 



(17) 



(18) 



Note that via Eq. (0) this Green's function confirms the form of Eq. (0) for t ^ 0, as well as the boundary condition 
Eq. (p^). Next, we consider the case of Langevin dynamics. 



dU 



(19) 



with the standard Gaussian white noise termtj {£,iit)^j{t')) — 2^kBTSi,jS{t — t'). From the equations of motion one 
finds that, because the integrated noise term is Gaussianly distributed. 



G{X -^X';t) ex exp 






J2 (xl -x^ + fi d^u tf 



i=l 



(20) 



At short times this again leads to Eq. (|^, where however now 



a{t) = - — , [Langevin dynamics, i — > 0] 



(21) 



In the case of Langevin dynamics Eq. (f^ may be applied to the calculation of the force-force time auto-correlation 
function. This leads to an equation that in principle allows a calculation of (AA'2(i)) directly form the double 
partition function: The mean-square displacement in time t is given by (sum over repeated indices) 



{AX\t)) - [ dt' [ dt" lx,{t')X,{t") 
Jo Jo ^ 



(22) 



Since the noise terms are uncorrelated at different times, Eqs. ( p^ and ( |2^ ) imply 



dt^ 



(Ax^it)) = 2(i:,(o)x,(i)) = 2^I^ {d,uio)d,u{t)) . 



(23) 



From Eq. (Jq) the force-force time auto-correlation function is rewritten as 



{dM0)dMt)) - i ^ '"''''' 



52 J D{a, b) 



a-e-''^^^)] k'e'''^^^') 



-a(X-X')^ 



By partial integrations one finds 



{d,u{o)d^u[t)) 



dXdX' 
52 J D{a,b) 



{X - X'f e-"('^-^')'-Mt^(^)+t^(^')] 



a2 91nD 
52 9a 



Thus, the equation for a{t) is from Eqs. (11), (|23|) and (Ea 



d2 ^ 2a2x 51^^ 



(24) 



(25) 



(26) 



The expansion of Eq. (E^) is straightforward, though tedious. 

In order to check the validity of Eq. (H) a simple system obeying Langevin dynamics was studied numerically. 
The system was chosen to be so simple that the integral of Eq. (pt) may be evaluated " exactly" , thus avoiding the 
noise of Monte Carlo simulations. No attempts were made to verify that Eq. (^6|) gives the correct a{t). Instead, 
the following procedure was followed. At a number of fixed a-values 5(a) was found from the requirement that the 
average joint potential energy is 2 ([/). Then the mean-square displacement was evaluated for each a from Eq. ( |ll| ) 
and also as function of time from the dynamical simulations, allowing an identification of the times corresponding 
to the fixed a-values. Finally, the time auto-correlation function was calculated from Eq. (pf) at the fixed a-values. 
Figure 2 shows the results for the foUoimng time auto-correlation function (^X^{0)X^{t)^ for a particle in a double well 
potential obeying LangevirL-dynamicsEd. The full curve is the exact time auto-correlation function found by solving 
the Smoluchowski equationE3 and the dots give the prediction of Eq. (pi) . Results are shown for the cases (3 = 2 and 
/3 = 8 in dimensionless units. There is good agreement between prediction and simulation. 

To conclude, in this paper a statistical mechanical approximation for the calculation of time auto-correlation 
functions has been proposed. The formalism assumes a knowledge of the mean-square displacement in configuration 
space as a function of time, the mean-square displacement acting as a "molecular clock". The remaining "spatial" 
auto-correlation is a "double" canonical average (Eq. (H)). Note that the "double partition function" , D(a, 6), contains 
the ordinary partition function (of configuration space) as a special case: Z^{f3) = £'(0, /?). 

The approximation proposed is only useful if (^AX'^{t)) is known. Experimentally, this quantity is accessible via 
the intermediate incoherent scattering function. In some cases phenomenological estimates of the mean-square dis- 
placement as function of time may be given. Thus, for hopping in a rugged energy landscape where all minima are 
equal, the mean-square displacement is universal at low temperatures iexcept for trivial scalings), i.e., indepen- 
dent of the barrier height probability distribution. This has been shownlij by the effective medium approximation 
and by comraiter simulations of the frequency-dependent conductivity (which is simply related to the mean-square 
displacementO) . Finally, there is the possibility that the mean-square displacement majj be found approximately 
from Eq. ([26^ if Langevin dynamics is assumed, as is common, e.g., in polymer dynamicsEJ. 

Equation (|26|) signals that Langevin dynamics plays a special role in the proposed scheme for calculation of time 
auto-correlation functions. A question of considerably interest is if ap_,when Langevin dynamics can be expected to 
give the same time auto-correlation functions as Newtonian dynamicsEj. If the ansatz proposed above is correct, two 
different dynamics give the same time auto-correlation functions for any quantity defined in configuration space, if 
the two dynamics just give the same mean-square displacement as function of time. In this way the ansatz provides 
a mechanism for the consistency of any two types of dynamics. 
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There are two figures. Figures and figure caption are available from author by snail mail, please write an E-mail to 
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